Search for new physics with neutrinos at Radioactive Ion Beam facihties 
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We propose applications of Radioactive Ion Beam facilities to investigate physics beyond the 
Standard Model. In particular, we focus on the possible measurement of coherent neutrino-nucleus 
scattering and on a search for sterile neutrinos, by means of a low energy beta-beam with a Lorentz 
boost factor 7 ~ 1. In the considered setup the collected radioactive ions are sent inside a 47r 
detector. For the first application we provide the number of events associated with neutrino-nucleus 
coherent scattering, when the detector is filled in with a noble liquid. For the sterile search we 
consider that the spherical detector is filled in with a liquid scintillator, and that the neutrino 
detection channel is inverse-beta decay. We provide the exclusion curves for the sterile neutrino 
mixing parameters, based upon the 3-1-1 formalism, depending upon the achievable ion intensity. 
Our results are obtained both from total rates, and including spectral information with binning in 
energy and in distance. The proposed experiment represents a possible alternative to clarify the 
current anomalies observed in neutrino experiments. 
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I. INTRODUCTION 

A wealth of experimental results on neutrino oseillations 
have been gathered sinee the neutrino oseillation discov- 
ery. Currently most of the data from aecelerators, reac- 
tors and the sun are interpreted within the theoretical 
framework of three active neutrino flavors involving the 
Maki-Nakagawa-Sakata-Pontecorvo matrix relating the 
flavor to the mass basis. In this case the number of un- 
known parameters is limited to three angles and three 
(one Dirac and two Majorana) phases, most of which has 
been determined. Recently, the T2K collaboration has 
found an indication for a non-zero value for the still un- 
known neutrino mixing angle ^13, at 2.5cr [1]. A non-zero 
6'i3 is also consistent with the first Double-Chooz results 
[2]. New results on the third neutrino mixing angle have 
recently been obtained by the Daya-Bay [3] and RENO 
[4] collaborations. The most precise measurement is cur- 
rently sin^ 26*13 == 0.092 ± 0.016(stot) ± 0.005(s?/sf), at 
5.2 a from Daya-Bay [3]. Note that a combined analysis 
had previously favoured a non-zero 6*13 (see e.g. [5, 6]). 
Beyond the intrinsic theoretical interest of knowing the 
last mixing angle value, its determination is a key step 
for setting up a strategy to search for leptonic CP viola- 
tion. With an upgrade of the T2K and NOi^A accelerator 
experiments, a (small) fraction of the Dirac 5 values can 
be explored [7] . The coverage of most of the Dirac phase 
values can be attained only with next generation experi- 
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ments including superbeams or beta-beams (see e.g. [8]). 
The DAEDALUS project constitutes an interesting alter- 
native [9] . Other open questions concerning fundamental 
neutrino properties include the neutrino mass scale, for 
which the KATRIN experiment should deliver results in 
the coming years [10], the neutrino mass hierarchy, the 
Majorana versus Dirac nature of neutrinos, and the pos- 
sible existence of sterile neutrinos. 

Besides the essential information gathered from ter- 
restrial experiments, neutrino properties have an impor- 
tant impact on astrophysical and cosmological observa- 
tions. Numerous examples exist in the literature showing 
that information can be extracted on unknown neutrino 
properties, or discussing their implications on a variety 
of phenomena like for example (stellar and primordial) 
nucleosynthesis processes. Recently it has been shown 
e.g. that a non-zero CP violating Dirac phase might 
have an impact in core-collapse supernovae [11, 12], or 
on Big-Bang nucleosynthesis [13]. Numerous studies have 
investigated the effects of sterile neutrinos, e.g. on the 
r-process (such as Ref. [14]) or on the primordial light el- 
ement abundances, like in Ref. [15]. Recent cosmological 
constraints on sterile neutrinos can be found in [16]. 

While neutrino oscillations are nowadays an estab- 
lished fact, several anomalies have recently been ob- 
served, that cannot be explained within the standard 
three active neutrino framework. First the MiniBooNE 
anti-neutrino and neutrino oscillation results are not fully 
understood, while an increased statistics should help to 
elucidate the low energy excess and the oscillation hy- 
pothesis [17]. This experiment, which was supposed to 
confirm/rule out LSND, has found an indication for neu- 
trino oscillations at a Am^ of about 1 eV^ both in the 



antineutrino channel, using dccay-at-rcst muons [18], and 
the neutrino channel, based upon decay-in-flight pions 
[19]. Note that the KARMEN experiment employing a 
similar neutrino source and detector has found no indica- 
tion for oscillations and has excluded most of the LSND 
oscillation parameter region [20]. The second anomaly 
is known as the reactor "anomaly" [21]. Indeed a recent 
reevaluation of the electron anti-neutrino flux from re- 
actors has shown a shift in the flux renormalization by 
3% [22], compared to the previous predictions. The re- 
analysis of the reactor experiments, using this new flux, 
has shown a significant inconsistency with the three neu- 
trino oscillation hypothesis. Finally, some years ago, 
the GALLEX and SAGE experiments pointed out an 
anomaly in the neutrino flux measured by putting an 
intense static ^^^Ar and ^^Cr sources inside their detec- 
tors. This is referred to as the Gallium anomaly. Ref.[23] 
has performed a detailed analysis including the 5-10% 
uncertainty on the corresponding neutrino-nucleus cross 
sections, showing that the Gallium anomaly statistical 
significance is at the level of 3cr. 

Currently the ensemble of the accelerator, reactor and 
Gallium anomalies are the object of debate and have trig- 
gered an intense investigation. The possible interpreta- 
tions exploit for example one or more sterile neutrinos, 
such as in [24], a combination of sterile neutrinos with 
non-standard interactions like in Ref. [25], while none of 
the proposed explanations so far provides a comprehen- 
sive understanding of all the data. Numerous propos- 
als are being put forward to confirm/rule out possibili- 
ties [26-30]. Among these, Ref. [26] proposes to put in- 
tense radioactive sources inside neutrino detectors, while 
Ref. [31] has pointed out the possibility to use intense ion 
sources produced at nuclear facilities. It is clear that in- 
dependent and aimed experiments are necessary, to clar- 
ify the present situation. 

Ref. [32] has proposed the idea of establishing a low en- 
ergy beta-beam^ facility, to dispose of neutrino beams 
in the 100 MeV energy range, based upon the beta de- 
cay of radioactive ions, with 7 w 1 (7 being the Lorentz 
factor) or with an ion boost 7 of typically 2 to 7^. In 
the first case, the neutrino fluxes are those of ions that 
decay-at-rest; while, in the second case, beams of vari- 
able average energy are obtained through a boost of the 
ions. The advantage of having such a facility is to dis- 
pose of pure (in flavor) and well known electron neu- 
trino (or anti-neutrino) fluxes. The physical applications 



^ The beta-beam concept was first proposed by Zucchelli to estab- 
lish a facility for the search of Icptonic CP violation [33]. For 
a discussion on the status of the feasibility of beta-beam facili- 
ties see e.g. [34]. Note that Ref. [35] has proposed a method to 
reach high Q-value ion intensities e.g. for *B and *Li, which is 
currently being investigated. 

■^ Obviously even larger ion boosts, around 10 or more, would be 
of interest. The numbers of 2-7 quoted in the available literature 
on low energy beta-beams was figured out to keep, in particular, 
the storage ring of small size. 



cover neutrino-nucleus interaction and fundamental in- 
teractions studies, oscillation searches and core-collapse 
supernova physics, as pointed out in [32] (for a review see 
[8] ) . These issues have been investigated in detail in a se- 
ries of works, including neutrino-nucleus interaction as- 
pects in [36-39] , a measurement of the neutrino magnetic 
moment in [40], a test of the Conserved- Vector-Current 
hypothesis [41], a measurement of the Weinberg angle at 
low momentum transfer in [42] , to non-standard interac- 
tions in [43] and of coherent neutrino-nucleus scattering 
in [44] , the oscillation towards sterile in [45] , an interpo- 
lation method to extract the neutrino fluxes from super- 
novae in [46]. In [47] the connection between neutrino- 
nucleus interaction and double-beta decay is pointed out 
in relation with a low energy beta-beam. IVIost of these 
applications are based on stored boosted ions. In Ref. [40] 
we have considered the configuration with 7^1, with 
radioactive ions sent to a target inside a Att detector for 
the search of the neutrino magnetic moment. Note that 
Ref. [31] has taken the same configuration for a sterile 
neutrino search. 

In this paper we consider a low energy beta-beam with 
7 ?a 1 . We consider that the ions are injected into a target 
inside a 47r detector. The purpose is to use the resulting 
very low energy neutrino flux to search for new physics. 
Here we explore two applications: the search for ster- 
ile neutrinos and a coherent neutrino-nucleus scattering 
measurement. We show that, depending on the ion inten- 
sity, a coherent neutrino-nucleus scattering measurement 
could be performed. Then, we focus on the search for 
one sterile neutrino in a 3-|-l-neutrino flavor framework 
and present exclusion plots for the sterile neutrino mix- 
ing parameters. The manuscript is structured as follows. 
We present our framework in Section II, while our nu- 
merical results are given in Section III. Section IV is a 
conclusion. 



II. GENERAL FRAMEWORK 

A. Possible setups and corresponding neutrino 
fluxes 

Radioactive ion beam facilities produce intense ra- 
dioactive ions decaying through beta-decay, or electron 
capture. Since specific radioactive ions can be selected, 
a pure electron (anti)neutrino flux can be obtained. As 
first proposed in [32] , there exists two alternative ways to 
produce low energy neutrinos (Figure 1). In the first sce- 
nario the decaying ions are stored inside a storage ring, 
while the emitted (anti)neutrinos are detected in a de- 
tector located close to the storage ring. If the stored 
ions are boosted, the corresponding neutrino spectra 
have variable energy with the average energy given by 
{E^) w 7(5/2, with Q being the Q-value of the beta- 
decaying nucleus. Depending on the application envis- 
aged, the neutrino fluxes can be tuned by appropriately 
choosing the Lorentz boost and a high/low Q-value ion. 
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The neutrino event rate detected on a spherical surface 
segment of width dr, located at a distance r from the 
center of the detector, is given by 



FIG. 1. (Color online) Two possible scenarios to produce low 
energy neutrinos from radioactive ions at Lorentz boost of 1. 
Left figure: The radioactive ions decay while circulating in a 
storage ring. The neutrino detector is installed close to the 
ring. Right figure: The radioactive ions are injected into a 
cavity inside a spherical detector. 



In the case the ions are not boosted (7 « 1), storing the 
ions in a small storage ring is a possibility as well. An 
example is furnished by the storage ring facility currently 
proposed at HIE-ISOLDE at CERN [48]. While for this 
specific storage ring the number of stored ions is limited, 
one can imagine the establishment of a small ring at one 
of the future intense radioactive ion beam facilities, such 
as European EURISOL [49], or the US Facility for Rare 
Isotope Beams (FRIB). 

The second scenario to produce low energy neutrinos, 
consists in injecting the ions into a target placed inside 
the detector. It turns out that, as long as radioactive 
ions are slow (i.e. not accelerated to Lorentz boost val- 
ues above 1), such a scenario is much more efficient than 
the storage ring case. This is due to a geometrical effect 
since, only part of the produced (anti)neutrino fiux - and 
not the total flux - traverses the detector if the ions are 
stored in a storage ring. The average neutrino fiux at 
the detector is further reduced, compared to the injec- 
tion inside the detector case, if the detector cannot be 
located very close to the storage ring due to the back- 
ground shieldings and other necessary instrumentation. 

In the rest frame, the beta-decay of a nucleus produces 
the following (anti-)neutrino fiux as a function of neu- 
trino energy 



Af{E,) = r^El E, ^JEI - ml F{Z, E,) Q{E, - m,) , 

(1) 
where / can be deduced from the measured /ii /2-value. 
The quantities appearing in the above expression are the 
emitted Icpton (electron or positron) energy E^^ = Q — E^, 
and the Fermi function F{Z, E^), which accounts for the 
Coulomb modification of the spectrum, with Z being the 
ion nuclear charge. 

We consider that radioactive ions are produced and 
injected into a target with an intensity / of ions per sec- 
ond. This target is installed inside a cavity located at 
the center of a spherical detector. The corresponding 
(anti) neutrino fiux at the distance r from the target is 



^-m^^r)..iE.)^^.r-dr. (3) 



Here ai(Ey) is the neutrino capture cross section on the 
target material z, Mmoi is the average molar mass of the 
detector material, Ui the average number of nuclei of type 
i per one mole of the detector material, and p is its den- 
sity. 



B. A coherent neutrino- nucleus scattering 
measurement 



The measurement of coherent neutrino-nucleus scat- 
tering constitute a precision test of the Standard Model, 
including the possibility to probe the weak nuclear charge 
as well as various deviations from the SM predictions, due 
to new physics above the weak scale, or presence of the 
sterile neutrinos [28] . Coherent neutrino- nucleus scatter- 
ing is also important in the astrophysical context, e.g. 
for core-collapse supernova physics. 

Several proposals have been made to perform such a 
measurement, particularly nearby Spallation Sources [50, 
51]. Here we consider a setup as shown in Figure 1 right. 
The cross section for electron neutrino (or anti-neutrino) 
coherent scattering on a nucleus is given by [52, 53] 
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here Gp is the Fermi constant, M is the nuclear mass, T 
the nuclear recoil energy, F the ground state clastic form 
factor and 



QvK = iV- {1 - Asin^ 9w) Z 



(5) 



the weak nuclear charge, with A^ the number of neutrons, 
Z the number of protons and Oyy the weak mixing an- 
gle. For neutrino energies below 50 MeV the momentum 
transfer is small and the form factor is close to unity 
F ^ 1. For the measurement of neutrino- nucleus coher- 
ent scattering, different types of liquids are being dis- 
cussed (see e.g. [50]). Here we take liquid neon as an 
example. 



C. A 3+1 sterile neutrino oscillation experiment 

In the present work we consider a sterile neutrino 
search within the 3+1 framework with three active 
neutrinos and one additional sterile neutrino. Besides 
the usual parameters of the Maki-Nakagawa-Sakata- 
Pontecorvo matrix, in this case the oscillation formula de- 
pends upon the neutrino mixing angle 9new and Am^^^ , 
considered to be much larger than Am^j^ ~ 2.4 • 10^'^ 
eV^. Implementing more complex scenarios with extra 
neutrinos is straightforward. The electron neutrino sur- 
vival probability for Pee {E^ , r) is given by Ref. [54] 



PeeiE^,r) = 1 - CDS'* 9new siu^ (26*13) siu^ 
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where a baseline of L < 2 km and neutrino energies E^ > 
2 MeV are assumed. Eqs.(l-3) are used to determine 
the unoscillated number of events, while Eq.(3) has to be 
multiplied by the neutrino survival probability Pge Eq.(6) 
in order to determine the number of oscillated events. 



1. Statistical analysis and backgrounds 



is performed over the subset of fitted parameters A not 
held fixed during the fitting, as well as over the 'pull' ^. 
The issue of the background levels is an important one. 
First of all we would like to stress that for the spherical 
detector setup (Figure 1, right) we will not have size- 
able beam'^ or implantation related backgrounds, since 
the ions implanted on the target are essentially at rest^. 
From previous experience with reactor experiments the 
primary sources being environmental and geoneutrino 
backgrounds, which one could deal with in a relative sim- 
ple fashion by implementing an energy threshold around 
4 MeV in the neutrino energy. The real problem are at- 
mospheric muons since one gets backgrounds via their 
high energy neutrons or their spallation left overs (Li-9). 
In order to reduce this kind of background signal to a 
tolerable level one needs a detector design similar to the 
one of, for instance, Daya Bay [57], where in addition to 
a rock overburden of at least 98 m (or 260 mwe), a 20 
ton liquid target must be surrounded by a 20 ton gamma 
catcher and 40 ton of buffer volume. Under these condi- 
tions the noise levels can be reduced to 6 events per ton 
per year. Therefore, to achieve such a low background 
event rate for the measurement under consideration here, 
it is necessary to locate the experiment well underground 
and surround the detector with appropriate shielding. In 
our calculations we assume that this can be done, and 
take as a reference value for the background, 6 events 
per ton per year. 



We present sensitivity plots obtained with the follow- 
ing procedure. We deal with systematic uncertainties 
inherent to the experimental setup by implementing the 
systematics directly into the statistical analysis by the 
use of the pull approach (see for instance [55, 56]). Un- 
less otherwise stated we bin our data in energy as well as 
in R intervals of equal spacing, R being the distance from 
the center of the detector, the x^ function being defined 
as 



X = mm^ J 



E^ 




(7) 
where the sum runs over energy and R-bins. As is cus- 
tomary in this type of analysis a superscript t is used to 
denote the predicted number of events TV*^, while a su- 
perscript / is employed to denote the number of events 
obtained from the fitting to the simulated data, n{- . The 
systematic error is indicated by tt. It enters the analysis 
through the definition: 
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n*, x(l + 7rC). 



(8) 



The respective true and fitted number of events, n* and 
n^-, are functions of their corresponding oscillating pa- 
rameters, but nonetheless the marginalization in Eq. (7) 



III. NUMERICAL RESULTS 

To produce low energy neutrinos both /3+ and /3~ de- 
caying ions can be considered as electron neutrino and 
anti- neutrino emitters respectively^. The choice of the 
ions depends on the achievable intensities, the half-lives 
and Q-values. Obviously the half-lives should lie in an 
appropriate range between short and long to make ex- 
periments feasible, so that typically half-lives in the 1 s 
range seem to be a good choice. On the other hand, high 
Q-values help increasing the total number of events as 
well as improving the signal-to-background ratio. 



^ Note that for the storage ring configuration (Figure 1, left) some 
background might be induced by the daughter nuclei colliding 
with the storage ring. Such a background can in principle be 
suppressed by putting an appropriate shielding around the stor- 
age ring. 

* Note that the situation here is very different from the one where 
the ions are boosted at high 7. In this case there is again no 
beam related background (a known advantage of the beta-beam 
concept) while there is radioactivity induced in the storage ring 
arising, e.g. from collisions of the stable daughter nuclei on the 
ring. In this case a shielding is necessary to suppress the related 
backgrounds and has been considered in the previous literature 
on low energy beta-beams with 7 larget than 1. 

"" Note that electron-capture neutrino beams have been considered 
in [58]. 



TABLE I. Beta-decay properties of the ions considered in our 
proposal: r is the decay lifetime, B™^'' is the end-point energy. 



Ion 

2 He 

iLi 

iLi 

IB 

?B 



Decay Daughter (State) 

/3- iLi(l+,0) 

/3- tBe(2+,0) 

/3- lBe(0+,0) 

P+ !Be(2+,0) 

13+ lBe(0+,0) 



r(ms) 

806.7 

838 

838 

770 

770 



£"=''= (MeV) 

3.5078 

13.103 

16.003 

15.079 

17.979 



Table I shows the candidate ions that we have been 
considering here, as typical examples. Note that there ex- 
ist a number of other promising radioactive ions, such as 
e.g. ®He or -'^^N [59]. As far as ®Li and *B are concerned, 
they decay mainly into a broad ^Be J'^ = 2+ excited 
state, at 3.03 MeV, therefore having Q-values centered 
at 13.1 MeV and 15.1 MeV, respectively. Nevertheless, 
due to the broadness of the final state, the neutrino spec- 
trum is extended well above the energy associated to the 
centered Q-value. To evaluate qualitatively this effect, we 
present results for two decay modes ^ : i) 100% branch- 
ing ratio to the ^Be ground state; ii) 100% branching 
ratio to a narrow excited state at 3.03 MeV. In a real 
experiment the actual results will fall in-between these 
two "extreme" cases. 

For the ion intensity, we assume 10^"^ ions per second. 
Instead of taking this parameter as a tunable number 
(as sometimes done in the literature), here we consider 
values, that can in principle be achievable at next gen- 
eration radioactive ion beam facilities. The predictions 
we present are obtained by taking 1 year = 10^ s and a 
100% efficiency of the detectors. 



A. Expected i^-nucleus coherent scattering events 

For the coherent neutrino-nucleus scattering measure- 
ment both electron neutrino (/?+) and anti- neutrino (/?^) 
emitters can be used. We have considered ^Li and ®B and 
their two decay modes (Table I). We take a 1 ton spher- 
ical liquid-neon detector, where the ions are injected in- 
side a central cavity having 20 cm radius (Figure 1 right). 
While other target nuclei are obviously possible, liquid 
neon is taken as an example. We would like to emphasise 
that detailed background simulations have already been 
done e.g. for the CLEAR experiment proposed at the 
SNS spallation source facility [50]. The shielding envis- 
aged for such a detector has been shown to reduce back- 
grounds nearby spallation sources at a negligible level. 
We expect that a similar reduction can be reached by 
putting the detector underground and/or using an ap- 
propriate shielding. However reaching very low nuclear 



TABLE II. Coherent i^-nucleus scattering: expected number 
of events for the two candidate ions considered and the setup 
of Figure 1 right (here 1 year — 10^ s). The maximal neu- 
trino energy is denoted by i?™^", the nuclear recoil detection 
threshold by Tmin. The results correspond to an intensity of 
10^^ ions/s. 

Ion Decay Target £?'^''(MeV) Tmin(keV) Events/ton/year 
iLi 13- Ne 13.103 10 192 

iLi /?- 

8 



B /3^ 
iB P^ 



Ne 
Ne 
Ne 



16.003 
15.079 
17.979 



10 
10 
10 



1373 



3047 



Note that an accurate neutrino spectrum might be obtained by 
describing ^Be final state continuum, when taking into account 
the delayed-a spectrum [60]. 



recoils is challenging for background issues (see, for ex- 
ample Figs. 10 and 11, of Ref. [50]) and for technical 
features, like light quenching (for a discussion see e.g. 
Ref. [61]). Although very optimistic, a sensitivity thresh- 
old of 10 keV will be assumed. A higher energy threshold 
choice can be an option in the kind of proposals discussed 
here, only if much higher ion intensities can be attained. 
Table II presents the number of expected events associ- 
ated to electron (anti)-neutrino scattering on neon. One 
can see that, despite the fact that coherence enhances the 
cross-section relative to other type of processes, the low 
energy range covered by the neutrino flux in this work 
makes the number of events still rather small, compared 
to the ones attainable with, e.g. the Michel spectrum of 
decay-at-rest muons produced at spallation sources [50]. 
Clearly, the feasibility of the measurement we are propos- 
ing strongly depends upon the achievable ion intensities 
and the lowest measurable nuclear recoil in the detec- 
tor. Figure 2 presents our predictions for the number of 
events for the two candidate ions considered, as a func- 
tion of the minimum measurable nuclear recoil. Note that 
it is straightforward to scale our rates (Table II) with the 
ion intensity collected at the center of the 47r detector, 
or to take into account effects such as light quenching. 
One can see the strong sensitivity of the results to the 
maximal neutrino energy depending on the Q— value of 
the ions. 



B. Active-to-sterile neutrino oscillation exclusion 
plots 

For the active-to-sterile neutrino oscillation search, we 
shall consider two types of /3~ decaying ions (Tabic I). 
First, because of its very well known aspects relating to 
its production and management, it is worth taking a look 
at the physics reach of a facility based on *'He. Its low 
Q-value yields fewer number of counts and, as we will 
show, this hinders the potential of a setup exploiting this 
ion, instead of one based upon a high Q-value ion choice. 
Our proposal for the search of sterile neutrinos is mainly 
based on the properties of ^Li, for which we assume two 
extreme cases, indicated as ^Li-16MeV and *Li-13MeV 
'ions'. 
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FIG. 2. Coherent i^-nucleus scattering: expected number of 
events, as a function of the minimal nuclear recoil for the 
experimental setup of Figure 1 right. The *B and *Li ion 
intensity is fixed at 10^''ions/s (here 1 year — 10^ s). 



Our choice of main setup has been dictated by an anal- 
ysis of the performance of the two possible configura- 
tions shown in Figure 1. In both cases, the considered 
detector is filled in with a liquid scintillator^. The elec- 
tron anti-neutrino detection channel is inverse beta-decay 
Pg + p — !• n + e+. A good signal-to-background ratio 
can be obtained by the addition of Gadolinium and the 
subsequent detection of the 8 MeV prompt gamma-rays 
produced by the neutron capture. 

For the case of the active-sterile neutrino oscillation 
hypothesis, under consideration here, we have chosen to 
present the results of our simulations by means of ex- 
clusion plots based upon Eq.(7). The exclusion plots for 
the active-to-sterile oscillation parameter sin^(26'now) are 
obtained by additionally fixing sin^ (26*13)® to the best fit 
value of Ref. [62], namely sin^(26'i3) = 0.051. Recently 
the Daya-Bay collaboration has measured the third neu- 
trino mixing angle at 5.2 a to be sin^ 26'i3 = 0.092 ± 
0.016{stat) ± 0.005{syst) [3]. We have checked that the 
plots presented here show no appreciable changes if the 
Daya-Bay value is used. The plots show the oscillation 
parameter space region where our setup is expected to be 
sensitive to the detection of active-sterile neutrino oscilla- 
tions. In all our calculations the considered running time 
of the experiment is of 5 years. Unless contrarily stated, 
we fix the systematic error to tt = 1%^ in all the analy- 
sis presented hereafter; while we will show how our main 
results change if a larger systematic error is considered. 
We shall compare the sterile neutrino oscillation param- 



^ We take as an example CisHig with a density of p = 988 kg/m^. 

* Note that the exclusion curves change little if one fixes the third 
neutrino mixing angle to zero. 

^ Note that it is not our goal to discuss how low systematic er- 
rors can be achieved in the actual experiment. For example, in 
Ref. [63] the authors discuss how this can be done in a short- 
baseline experiment at a neutrino factory. 



eter regions that can be covered with our experimental 
setup, to the allowed regions presented in the analysis of 
Ref. [21], based on reactor neutrino experimental data 
cumulated so far. 

Figure 3 presents exclusion plots calculated from a sta- 
tistical analysis of the data using total rates. The facility 
is based on ®Li-13 MeV decaying ions. The aim of the fig- 
ure is to compare the results obtained for the two experi- 
mental setups of Figure 1. Note that for the specific case 
of the storage ring only, we assume an intensity of 10^^ 
ions/s, having in mind a facility like HIE-ISOLDE (al- 
though the stored ion intensity is expected to be smaller 
[48]). Such an intensity should be attainable in a storage 
ring nearby the EURISOL facility [49]. For the setup 
geometry, following the TSR proposal for HIE-ISOLDE, 
we take a square storage ring with straight sections of 
61.6 m length and a 1 kton cubic^° detector at the stor- 
age ring center. (Such a geometry leaves 3 meters space 
between the detector and the storage ring straight sec- 
tions [64, 65].) Note that, the number of expected events 
and, thus the exclusion plots, strongly depends on the 
setup geometry. For a large detector, as considered here, 
placing it in the center of the storage ring represents the 
optimal scenario (if such a detector is located along one 
storage ring straight section, the event number is reduced 
by almost a factor of 4) . As expected, although the de- 
tector is of 20 tons only, the performance obtained by 
sending the ions inside a 47r detector is superior to the 
storage ring one^^ with respect to the coverage of the 
shaded region identified by the reactor anomaly. On the 
other hand, the storage ring setup has a better sensitivity 
to small Arn^. 

In addition to the aforementioned geometric advantages, 
the spherical detector setup benefits from the fact that 
neutrino source is very close to the active detector mate- 
rial (we remind that the radius cavity is of 20 cm only) . 
From now on, all the results we present will correspond 
to the spherical detector setup. 

Figures 4 and 5 show exclusion plots constructed from 
total rates, and from a spectral (binned) analysis of the 
simulated data, respectively. For the binned case, 8 
energy-bins and 8 R-bins of equal size are considered. 
Note that, for the binned case, the corresponding back- 
ground in the detector is scaled accordingly (it grows 
radially as i?^). Results for three ion sources are shown: 
the *^Li-16 MeV, the *Li-13 McV and ^He cases. The low 
Q-value of the helium ions clearly hinders the sensitiv- 
ity of this setup, making it clearly inferior to the lithium 
ion source case. Notice the slight difference between the 
two ®Li cases, which is only marginally enhanced for the 
binned case for large A'ni^^^{>7 eV^). (Small) correc- 
tions from ions decaying to the ground state of the daugh- 



^^ The detector base has a size of 9.4 X 9.4 m and the height is of 
11.3 m. Half of the detector is located below the storage ring, 
and half above. 

^^ Note that this is also due to the higher ion intensity. 
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FIG. 3. Exclusion plots for the sterile neutrino mixing param- 
eters from an analysis of the data including only total rates. 
The results are obtained by considering a 3+1 neutrino oscil- 
lation formalism. The contours shown are for a C.L. of 99% (2 
d.o.f.). The two setups are those of Figure 1. The solid (red) 
line corresponds to the 4n detector surrounding the source; 
while the dashed-dotted (blue) line corresponds to the detec- 
tor place at the center of the square storage ring. The ion 
intensities are of 10^^ ions/s for the storage ring and of 10^"^ 
ions/s for the spherical detector (see text). The source is *Li 
ions decaying mainly to the first excited state of the daughter 
nucleus (maximal neutrino energy 13 MeV). In both cases a 
5 years running time is assumed. For comparison the shaded 
region represents the 99% C.L. inclusion domain, given by 
the combination of reactor neutrino experiments and other 
sources (adapted from Fig. 8 of Ref. [21]). 
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FIG. 4. Exclusion plots for the sterile neutrino mixing pa- 
rameters from an analysis of the data including only total 
rates. The contours shown are for a C.L. of 99% (2 d.o.f.). 
The results correspond to choosing different ions: a source 
of *Li ions decaying mainly to the first excited state of the 
daughter nucleus (solid red line), *Li ions decaying mainly 
to the ground state of the daughter nucleus (dashed-dotted 
blue line), or a source of ®He ions (dashed brown line). For 
comparison the shaded region represents the 99% C.L. inclu- 
sion domain, given by the combination of reactor neutrino 
experiments and other sources (adapted from Fig. 8 of Ref. 
[21]). 



FIG. 5. Exclusion plots for the sterile neutrino mixing pa- 
rameters, with binned analysis of the simulated data. Both 
binning in neutrino energy and in distance within the detec- 
tor. The contours shown are for a C.L. of 99% (2 d.o.f.). 
The results correspond to choosing different ions: a source 
of *Li ions decaying mainly to the first excited state of the 
daughter nucleus (solid red line), *Li ions decaying mainly 
to the ground state of the daughter nucleus (dashed-dotted 
blue line), or a source of ®He ions (dashed brown line). For 
comparison the shaded region represents the 99% C.L. inclu- 
sion domain, given by the combination of reactor neutrino 
experiments and other sources (adapted from Fig. 8 of Ref. 
[21]). 



ter nucleus are thus expected to be important only in the 
large Arn^^^ case. The results of Figure 5 show the im- 
portance of an appropriate binning. 

For comparison we have also included, in these figures, 
shaded regions corresponding to the 99% C.L. inclusion 
domains identified by the combination of data from the 
reactor neutrino experiments and other sources as de- 
scribed in, and here adapted from. Fig. 8 of Ref. [21]. One 
can see that the proposal investigated here would allow 
to cover most of the active-sterile oscillation parameter 
region. On the other hand we recall that the presented 
exclusion curves have the following simple physical mean- 
ing: an actual measurement lying inside the curve (to the 
upper-right of the curve) represents definite evidence in 
favour of the corresponding hypothesis, in our case, ac- 
tive neutrinos oscillating into sterile ones. In this man- 
ner, from Figure 5 one sees that the shaded region is out 
of reach if one uses 10^'^ ^He/s; whereas using 10^'^ ®Li/s 
one can cover around 70%-75% of the currently allowed 
region. 

We would like to discuss now the impact of the cho- 
sen ion intensities on the setup performance. Figure 6 
shows how the exclusion plots (and the coverage of the 
allowed region) changes when varying the ion intensity. 
In particular, the physics potential, relative to our main 
setup with lO"'^'^ ions/s, is seen to diminish (increase) by 
changing the intensity in one order of magnitude. This 
speaks of the high level of influence, that achieving good 
ion production levels nearby future radioactive ion beam 
facilities have, upon these type of experimental searches. 
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FIG. 6. Exclusion plots with binned analysis of the simu- 
lated data, obtained by varying the ion intensity; 10^* ions/s 
(dashed-dotted), 10^'^ ions/s (solid) and 10^^ ions/s (dashed). 
The red lines correspond to the source of *Li ions decaying 
mainly to the first excited state of the daughter nucleus (max- 
imal neutrino energy 13 MeV); while the brown lines corre- 
spond to the source of ®He ions. The contours shown are for 
a C.L. of 99% (2 d.o.f.). For comparison the shaded region 
represents the 99% C.L. inclusion domain, given by the com- 
bination of reactor neutrino experiments and other sources 
(adapted from Fig. 8 of Ref. [21]). 



FIG. 7. Exclusion plots with binned analysis of the simu- 
lated data, obtained for the 10^"^ ions/s with *Li ions decay- 
ing mainly to a 13 MeV excited state of the daughter nucleus. 
The figure shows the impact on the exclusion curves of dif- 
ferent levels of systematic errors, namely 1% (solid, red), 2% 
(dotted, black), 5% (dashed, brown) and 10% (dashed-dotted, 
blue). The contours shown are for a C.L. of 99% (2 d.o.f.). 
For comparison the shaded region represents the 99% C.L. in- 
clusion domain, given by the combination of reactor neutrino 
experiments and other sources (adapted from Fig. 8 of Ref. 
[21]). 



Finally, the sensitivity of the proposed experiment 
might depend upon the achieved systematic errors. To 
show their impact, we present exclusion curves based 
upon a binned analysis for sterile neutrino mixing pa- 
rameters, for different levels of systematic errors. Fig- 
ures 7 and 8 show the impact of 1%, 2%, 5% and 10% 
systematic error on the exclusion curves, for 10^'^ and 
10^^ ^Li/s, respectively. One can see the important im- 
pact that reaching low systematic errors has, especially 
for large Am^ and small mixing angle. 
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IV. CONCLUSIONS 

Future intense radioactive ion beam facilities can offer an 
unique opportunity to perform searches for beyond the 
Standard Model physics, using low energy neutrino fluxes 
from beta-decaying ions. Here we consider two configu- 
rations, where either the ions are stored in a storage ring, 
or they are sent into a target inside a spherical detector, 
filled in either with a noble liquid, or with a scintilla- 
tor (with the addition of Gadolinium). Our results show 
that, as long as the ions are not boosted, the spherical 
geometry scenario gives better results than the storage 
ring one. We have presented predictions for the expected 
events associated to a coherent neutrino-nucleus scatter- 
ing measurement. The realization of such an experiment 
heavily depends on the achievement of large ion intensi- 
ties and reaching challenging low energy nuclear recoils. 
The second option considered is a sterile neutrino search, 
that can be performed using electron anti-neutrino detec- 
tion through inverse beta-decay in a scintillator. We have 



FIG. 8. Exclusion plots with binned analysis of the simu- 
lated data, obtained for the 10^* ions/s with *Li ions decay- 
ing mainly to a 13 MeV excited state of the daughter nucleus. 
The figure shows the impact on the exclusion curves of dif- 
ferent levels of systematic errors, namely 1% (solid, red), 2% 
(dotted, black), 5% (dashed, brown) and 10% (dashed-dotted, 
blue). The contours shown are for a C.L. of 99% (2 d.o.f.). 
For comparison the shaded region represents the 99% C.L. in- 
clusion domain, given by the combination of reactor neutrino 
experiments and other sources (adapted from Fig. 8 of Ref. 
[21]). 



presented exclusion plots obtained from total rates and 
from analysis including spectral information (with bin- 
ning in neutrino energy and distance within the detector) 
of the simulated data. In particular, the binned analysis 
gives interesting results for ion intensities achievable at 
future radioactive ion beam facilities, like e.g. the EU- 
RISOL facility. Clearly the ion intensities achievable at 



such facilities arc lower than the MCi radioactive source 
considered in proposals like the one in Ref. [26] . However 
radioactive ion beam facilities offer the possibility to dis- 
pose of radioactive ions with different Q-values, allowing 
to cover different regions of the oscillation parameters. 
With our spherical setup, one can probe large squared- 
mass differences and rather small mixing angle values, 
associated with one sterile neutrino, in the 3+1 oscilla- 
tion framework. In particular, with the kind of setup 
we consider here, one could confirm/rule out the sterile 
neutrino hypothesis, as a possible explanation of the cur- 
rently debated reactor neutrino anomaly, and cover most 
of the corresponding parameter space region. 
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